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Abstract. We study multivariable Christoffel-Darboux kernels, which may be viewed as 
reproducing kernels for antisymmetric orthogonal polynomials, and also as correlation func- 
tions for products of characteristic polynomials of random Hermitian matrices. Using their 
' interpretation as reproducing kernels, we obtain simple proofs of Pfaffian and determinant 

formulas, as well as Schur polynomial expansions, for such kernels. In subsequent work, 
these results are applied in combinatorics (enumeration of marked shifted tableaux) and 
number theory (representation of integers as sums of squares). 
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• Dedicated to the memory of Vadim Kuznetsov 

1 Introduction 

The Christoffel-Darboux kernel plays an important role in the theory of one- variable orthogonal 
polynomials. In the present work, we study a multivariable extension, which can be viewed 
as a reproducing kernel for anti- symmetric polynomials. As is explained below, our original 
motivation came from a very special case, having applications in number theory (sums of squares) 
and combinatorics (tableaux enumeration) . More generally, this kind of kernels occur in random 
matrix theory as correlation functions for products of characteristic polynomials of random 
Hermitian matrices. The purpose of the present note is to highlight a number of useful identities 
for such kernels. Although, as we will make clear, the main results can be found in the literature, 
they are scattered in work belonging to different disciplines, so it seems useful to collect them 
in one place. Moreover, our proofs, with the interpretation as reproducing kernels, are new and 
conceptually very simple. 

We first recall something of the one-variable theory. Let 

/ i-> / f(x) d/j,(x) 
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be a linear functional defined on polynomials of one variable. We denote by V n the space of 
polynomials of degree at most n — 1. Assuming that the pairing 



(f,g) = J f(x)g(x)dfi(x) 

is non-degenerate on each V n (for instance, if it is positive definite), there exists a corresponding 
system (pk(x))kLo of monic orthogonal polynomials. We may then introduce the Christoffel- 
Darboux kernel 

which is the reproducing kernel of V n , that is, the unique function such that (y i— ► lif (a;, y)) £ 
and 

f(x) = J f(y)K(x,y)d f i(y), f E V n . 

The Christoffel-Darboux formula states that 

1 PnO)Pn-l(y) -Pn-l(^)Pn(?/) 



K(x,y) 



{Pn-l,Pn-l) X-y 



More generally, let V™, < m < n, denote the mth exterior power of V n . It will be identified 
with the space of antisymmetric polynomials f(x) = f(x±, . . . ,x m ) that are of degree at most 
n — 1 in each X{. Writing 

f(x)dn m (x) = — / f(xi,...,x m )d(i(xx)---dfj,(x m ), 
ml J 

we equip V™ with the pairing 

(f,g)v™ = J f(x)g(x)dfi m (x). 

Equivalently, in terms of the spanning vectors det (fj{xi)), fj 6 V n , 

l<i,j<m 

(det {fj{xi)) det (gj(xi))) = det ({fi,gj) Vn )- (1) 

Every element of V^ 71 is divisible by the polynomial 
A(x) = (xj-Xi). 

l<i<j<m 

The map / ^ //A is an isometry from V™ to the space W™, consisting of symmetric polyno- 
mials in xi, . . . , x m that are of degree at most n — mm. each Xi, equipped with the pairing 

(f,g)w™ = —. I f(x)g(x)A(x) 2 dfi(xx)---d(i(x m ). 
ml J 

We remark that, normalizing 

A(x) 2 d/i(xi) • • • d(JL{x m ) (2) 
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to a probability distribution, it defines an orthogonal polynomial ensemble (the term is sometimes 
used also for the more general weights |A(x)|^). Such ensembles are important in a variety of 
contexts, including the theory of random Hermitian matrices |12j . 
An orthogonal basis of V™ is given by (ps(%))sc\n], \s\=mi where 

Ps(x) = det (pj-iixi)). 

l<i<m, j£b 

Here and throughout, we write [n] = {1, . . . ,n}, and we assume that the columns are ordered 
in the natural way. Indeed, it follows from ([T|) that 

(ps,Pt) = S S t Y[{pi-l,Pi-l). (3) 

We denote by A(x)A(y)K m (x, y) the reproducing kernel of V™, that is, the unique element 
of y n m ® 1/ m such that 

f(x) = J f(y)A(x)A(y)K m (x,y)d f i m (y), f G V™. (4) 

Equivalently, K m (x,y) is the reproducing kernel of W™. 
It is easy to see that 

Indeed, it follows from (JTJ) and ([3|) that both sides satisfy (jl]). The equality of the two expressions 
can also be derived using the Cauchy-Binet formula. 

The following elegant integral formula was recently obtained by Strahov and Fyodorov [29]. 
In Section [2] we will give a simple proof, using the interpretation of the left-hand side as a 
reproducing kernel. 

Proposition 1 (Strahov and Fyodorov). One has 

K m (x,y) = — I Y[ (zj - w k ) A(w) 2 dfi n _ m (w), (6) 

n(pj_l,Pi_!) J l<J<2m 
(=1 l<fc<n— m 

where 

(z 1 ,...,z 2 m) = (x 1 ,...,x m ,y 1 ,...,y m ). (7) 

As a consequence, the obvious S m x S m x Z2 symmetry of K m (x,y) extends to a non-trivial 
52m symmetry. 

Corollary 1. The polynomial K m (z) = K m (zi, . . . , z<im) is symmetric in all its variables. 

The motivation for the work of Strahov and Fyodorov is an interpretation of the integral ([6|) 
as a correlation function for the product of characteristic polynomials of random Hermitian 
matrices. The case x = y is of particular interest, since Q may then be written as 



A(xi , . . . , Xm) K m (x\ , . . . , x mi x\ 
1 

n 

(n - m)\ n (Pi-i,Pi-i) 

i=l 



A(xi, . . . ,x n ) 2 dn(x m+ x) ■ ■ ■ dfi(x n ), (8) 
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which exhibits A(x) 2 K m (x, x) as a correlation function for the measure ([2]). In that case, the 
determinant formula (0) is classical, see Chapter 5]. In our opinion, our proof of Proposition!]] 
is more illuminating than the inductive proof usually given for the special case x = y. For 
applications of Proposition Q] and for further related results, see [HEJH]. 

An alternative determinant formula for the integral ([6]) may be obtained by combining two 
classical results of Christoffel [H Theorem 2.7.1] and Heine [9j Theorem 2.1.2], see also [4]. We 
will obtain it as a by-product of the proof of Proposition [TJ 

Proposition 2. In the notation above, 

1 ^ de X (Pn-m+j-lOi)) 

lZj "n/ v A(z) 

[[{Pn-i,Pn-i) 
i=l 

More generally, such a determinant formula holds for the integrals 



/ II i z j - w k) &{wf dii n {w), 



l<jr'<m 
l<fc<n 

but we focus on the case when m is even. 

The fact that the "two-point" determinant in ([5]) and the "one-point" determinant in Propo- 
sition [2] agree can also be derived from the work of Lascoux |15| Propositions 8.4.1 and 8.4.3], 
see |151 Exercise 8.33] for the integral formula in this context. Note that the reproducing pro- 
perty mentioned by Lascoux, and given explicitly in [171 Proposition 3], is of a different nature 
from @, pertaining to integration against one-variable polynomials. 

The special case of Proposition [2] obtained by subtracting the fcth row from the (m + k)th, 
for 1 < k < m, and then letting z m+ k — > Z/., gives the following formula for the correlation 
function (jSJ). 

Corollary 2. In the notation above, 

n (p„_„ ft ,_ i )A( I )' ls - jS2m V K- m+ i-iM. m + l<i<2m) 

1=1 

Next, we give Pfaffian formulas for K m . As is explained below, they can be deduced from 
results of Ishikawa and Wakayama [8], Lascoux [16j . and Okada [22]. Nevertheless, we will 
give an independent proof, using Corollary [U Recall that the Pfaffian of a skew-symmetric 
even-dimensional matrix is given by 



K^L (oy) = 2^! £ m*) ik* 

— ' J — 0"6i>2m *=1 

Proposition 3. For any choice of square roots 



i ■ 



K m (z) = ^ r ~— =r pfaff ((^/z~ - y/Ti)K(zi, zj)) . 

11 W z i ~ y Zi > l < i J< 2m 

l<i<j<2m 

Moreover, for any choice of Q such that 

d + cr 1 = * + 2 > (9) 

2m 

n cr -1 

g OT ( Z )= * =1 — pfaff ((Ci-Ci)^,^-))- 

11 ISj ~ SiJ l<«,j<2m 
l<i<j'<2m 
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Note that © implies 

zj -zt = (Cj - C0(i - CiCi)- (io) 

In the special case 2 = (a;, a;), choosing 

• • • ) y/z2m) = {~ \J x l j ■ ■ ■ j ~^/Xmi \J x l, • • • > V a^mjj 
(Cl) • • • >C2m) = (£f\ • • • j^ifli • • • ,£m), 

Proposition [3] reduces to special cases of the identity 



pfaff 



A B 
-B -A 



\ = (-l)h m ( m ~ l ) det{B - A). 



The general case seems to lie deeper. 

Proposition[3]is actually equivalent to PropositionHJbelow, which can be deduced from known 
results. Indeed, rewriting (jllap using [22], Theorem 4.2] and (jllb[) using the case n = m of [22l 
Theorem 4.7], we can easily see the resulting expressions to agree. More explicitly, this identity 
appears in [16] , with a simple proof. The Pfaf f ian (jllcp can be treated similarly, or else shown to 
agree with (fTTbl) by means of a result of Ishikawa and Wakayama Theorem 5.1], see Remark [3] 
below. 

We remark that the relevant results of [8], [16] and [22] are closely related to Sundquist's 
identities [30], see also [6] and [7]. Moreover, (lllaj) is related to the Izergin-Korepin determinant 
for the partition function of the six- vertex model [10], which, as well as the Pfaffians in (jlip . 
has applications to alternating sign matrices [131 [T4"l [23| [32] . 

Proposition 4. Let en and Zi, 1 < % < 2m, 6e /ree variables, and let Q be as in ([9]). Moreover, 
let S C [2m] 6e an arbitrary subset of cardinality m. Then, 

i m(m+ i)+£s / a a A 
-1 ses (zj-zi) det M 11a 

= ]I (v^+v^) Pfeff f ^~ V ) (nb) 

l<i</<2m l<i,i<2m + V*iJ 

2m 



i=l l<i<j<2m 



Note that Proposition [3] is a special case of Proposition 0] when = p n (zi) / p n ~i(zi) , 
{x\, . . . , x m } = {zi}i£s an d {yi, • • • , y m } = { z i}igs- Conversely, it is not hard to deduce Propo- 
sition H] from Proposition [3] using an interpolation argument. 

Remark 1. Strahov and Fyodorov found one-point and two-point determinant formulas for 
more general correlation functions than those of Proposition [TJ In the case of 



l<k<n 



their two-point formula is of the type (|lla|) . see |29t Proposition 4.2]. Thus, Proposition 0] 
implies Pfaffian formulas for this correlation function. 
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In [27], we need an elementary result on the expansion of the kernel K m (x,y) into Schur 
polynomials s\(x)s fJi (y). Since we have not found a suitable reference, we include it here. The 
proof is given in Section [H 

Proposition 5. One has 

m 

v-^ (-1V=1 

K m (x,y) = 2^ . cr nc e lr n J c *+i-2) s A(x)s / ,(?/), 

0<Am<---<Ai<n-m J J , 
0</i m < - <A»i<n~m j=l 



S = {Afc + 771 + 1 — /c; 1 < k < m}, T = {fi^ + m + 1 — k; 1 < k < ?n}, 

and 

Cfc = y x fc dfj,(x). 

Finally, let us describe our original motivation, which comes from applications that seem 
completely unrelated to random matrix theory. Motivated by the theory of affine superalgebras, 
Kac and Wakimoto [11] conjectured certain new formulas for the number of representations of an 
integer as the sum of 4m 2 or 4m(m + 1) triangular numbers. These conjectures were first proved 
by Milne [HI [201 EI] , and later by Zagier [31] . In [25] , we re-derived and generalized the Kac- 
Wakimoto identities using elliptic Pfaffian evaluations. Extension of this analysis from triangles 
to squares leads to formulas involving Schur Q -polynomials [18] evaluated at the point (1, . . . , 1). 
(More precisely, these polynomials are normally labelled by positive integer partitions. Here, 
we need an extension to the case when some indices are negative.) Later, we realized that the 
resulting sums of squares formulas are equivalent to those of Milne [21] . Seeing this is far from 
obvious and requires an identification of the Schur Q-polynomials with kernels K m (z), where the 
underlying polynomials Pk(%) are continuous dual Hahn polynomials. The key fact for obtaining 
this identification is the second part of Proposition [3j We refer to [26] for applications of the 
results above to Schur Q-polynomials and marked shifted tableaux, and to [27] for the relation 
to sums of squares. 



2 Proof of Propositions [T] and [2] 

Lemma 1. Let cf> : V™ — > V^~ m be defined by 
(<f>f)(x) = J f(y)A(y,x)d f i m (y). 

Then, 

|m(m+l)+ J2 s f-f 
(#s)(z) = (-1) seS [[(Pi-l,Pi-l) -Ps-(x), 

where S° = {1, . . . , n} \ S. 

A corresponding statement holds when V n is a general n-dimensional vector space and A an 
element of the one-dimensional space {V*) n [3j § 8.5], the map <f> : V™ — > (V*) n ~ m often being 
called Hodge star or Poincare isomorphism. For completeness, we provide a proof in the present 
setting. 



Multivariable Christoffel-Darboux Kernels 



7 



Proof. The Vandermonde determinant evaluation 

A(x) = del -.(pj-iixi)) (12) 

±<i,j<n 



gives 

„ m n 

r(i)-l(»i)> 



/m n 



<r: [m]—*S i=l i=m+l 

r: [n]->[n] 

where the sum is over bijections. By orthogonality, we may assume 

(r(l), . . . , r(n)) = (cr(l), . . . , er(m), p(l), ...,p(n- m)), 
with p a bijection [n — m] — > S c . One easily checks that 

|m(m+l)+ E s 
sgn(r) = (-1) »es sgn (cr)sgn(p), 

which gives indeed 

i m(m+ i )+ £ s i ™ 

= (-1) S6S 2^ ll(2V(i)-l>2Mi)-l> 

' cr: [m]-»Si=l 
n— m 

X J] Sgll (^ II Pp(i)~l( X i) 
p:[n~m]^S c i=l 

|m(m+l)+ £ « 
= (-1) ses 1 1 (pi_i s Pi_i> p^e (a;) . 

ies 



By iteration, it follows from Lemma [T] that 

n 

4, O 4> = (_!)-(—) Hip^p^) . idyrn . 



1=1 



By means of A(w,x) = (— l) m ( n m ) A(a;, it;), x G R"\ m , this fact can be expressed as 

f( x ) = H " I f(y)A(x,w)A(y,w) dfi m (y)dfi n - m (w), f G V™. 

n (Pi-i,Pi-i) 

i=l 

By the uniqueness of the reproducing kernel, it follows that 

A(x)A(y)K m (x,y) = — / A(x,w)A(y,w)dfi n - m (w). (13) 

il {Pi-i,Pi-l) 

i=l 

This is equivalent to Proposition [TJ 



Remark 2. The equation (|13j) can also be obtained as the special case I = n of the contraction 
formula 

A(x)A(y)K m (x,y) = ± — J ' y / A(x,w)A(y,w)K l (x,w,y,w) dni_ m (w), (14) 



(n — m)! 

< m < I < n. Conversely, (|14|) follows easily from (|13l) . 
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To prove Proposition [2l we note that 

A / \ A(x,y)A(x,w)A(y,w) 
A ^ W ^= A( X )A(y)A( W ) 

Applying this to (fT3j) gives 



1 1 

K m (x,y) = — -r-, r / A(x,y,w)A(w)dfi n _ m (w) 

n<w-i,ft-i> A(!B ' v) 

1 {4>A)(z) 



" TT/ \ A ( Z ) ' 

lI(Pi-l,P»-l) 

( = 1 

where, as in (J7|), z = (x,y). Next we observe that, by (fT2|) . A(w) = ps(w), with 5 = [n — m]. 
Lemma [J then gives indeed 

1 ^ de )o (Pn-m+j-lOi)) 

^m(x,y) = . 

n(Pn-i,Pn-i) 
i=l 

3 Proof of Proposition [3] 

Our main tool is the following elementary property of Pfaffians, which we learned from an 
unpublished manuscript of Eric Rains [23] . 

Lemma 2 (Rains). For arbitrary (aij)i<ij<2m> 

pfaff (a ij -a ji )= V (-1)*^ det (ay), 

where x(S) denotes the number of even elements in S. 

For completeness, we sketch Rains' proof. 
Proof. The left-hand side is given by 

^ m .. m 

^ sgn(ff)JJ(a (T(2i _ lM 2i ) -a (T(:hV( 2 i _ 1) ) = ^ ^ sgn(cT) JJ a CT(2 ;_ 1)|(7(2i) . 

Write 

m 

^ a o-(2i-l),(r(2i) = H a j,r(j), 
i=l ieS 

with S 1 = {cr(l), a(3), . . . ,o~(2m — 1)} and r a bijection S — > S c . Identifying r as an element 
of S m , using the natural orderings on S and S c , it is easy to check that sgn(<r) = (-l)*( s ) sgn(r). 
Since the map a t— > (5, r) is m! to one, we obtain indeed 

£ (-l)x(s)^ sgn(r) TJ G . r(i) . ■ 

SC[2m],|5|=m r igS 

The following identity appeared as [H Theorem A.l]. 
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Corollary 3 (Ishikawa and Wakayama). One has 

Zj Z{ 



pfaff 

l<i,j<2m \ a + °{ x i + x j) + cx i x j 
= {b 2 -ac)^™" 1 ) 



a + b(xi + xj) + cXiXj 



l<i<j<2m 

x ^2 (-l) x(S,) JJ Zj (xj - Xi)(a + b(xi + Xj) + CXiXj). 

SC[2m],\S\=m j<£S l<i<j<2m 

i,jeS or i,j£S 

This follows immediately from Lemma O by use of the Cauchy determinant 

n ( x j -*i)(vj -I/O 

2 _ nf ,\lm(m-l) l<i<3<™ 



det T7 x ) =(b ~ ac) 



l<y<m V« + /M.r ; + //,) + cXiUj J » (<j + ^ + ^ + ^ ' 

which is reduced to its more well-known special case a = 1, 6 = 0, c = —1 through the elementary 
identity 

a + b(xi + Xj) + CXiXj = - (ac — b 2 + (cSj + b)(cxj + 6)) . (15) 
The proof of Corollary [3] given in [8] is more complicated. 

Remark 3. The equality of (jllbp and (|llc[) can be deduced by applying Corollary [3] to both 
Pfaffians, using also IjlOlh 



We only need Corollary [3] in the case when X{ = Z{. Then, the Pfaff ian is given by 

(&2 _ ac )m(m-l) TT . 

i< l ij<2 m a + b ^ + x ^ + cx ^ 



Indeed, one may use (|15p to reduce oneself to the case a = 1, b = 0, c = —1, which is the 
Pfaff ian evaluation in [281 Proposition 2.3]. 

Corollary 4. One has 

^ (-l) x(s,) (xj - Xj)(a + 6(xi + Xj) + caciXj) 

SC[2m], |S|=m j^S l<i<j<2rn 

i,j£S or ij'^S 

= (6 2_ Gc) | m (m-l) JJ fo-^). 

l<i<jf<2m 

We are now ready to prove Proposition [3l By Lemma [21 we have in general 

pfaff {(a j -a i )K{z i ,Zj))= V (-l)x(^) TT det (K(z h Zj)), 
^^ 2m sQ2^\s\= m Jfs ieS ' j * S 

which, by Corollary [H equals 

K m ( Z ) yi (-i) x(s) n^- n (**-*)■ ( ie ) 

SC[2m],\S\=m j£S l<i<j<2m 

i,j£S or i,j£S 
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Consider first the case en = y/z^ Then, by the case a = c = 0, 6 = 1, X{ = ^fzi of Corollary [H 
the sum in (|16|) equals 

n - • 

l<i<j<2m 

This yields the first part of Proposition [3l Similarly, letting ai = Q and using (|10p , we can 
compute the sum in (|TB|) by the case a = 1, b = 0, c = —1, = £j of Corollary H] as 

2m 

nc ? m n &-<*)■ 

t=l l<i<j<2m 

This completes the proof of Proposition 

4 Proof of Proposition [5] 

When (efc)^ =1 is a basis of V n , let (es)gc[n], \S\=m be the corresponding basis of V™ defined by 

es(x) = det (ej(xi)). 
l<i<m,jeS J 

We then have the following general expansion formula. If ei(x) = fi(x) = pi-\(x), this is ([5]). 
Lemma 3. Let (efe)^ =1 and (fk)k=i ^ e arbitrary bases ofV n . Then, 

A(x)A(y)Jf m (x,»)= (-!) seS f£T /w/ m ^WHlz)- (17) 

S,TC[n] l<iJ<n Uei ' /j ' ;J 
|5|=|T|=m 

Proof. It is enough to show that 

e u(x) = eu(y)R(x,y)dfi m (y), U C [n], |J7| = m, 



where i? denotes the right-hand side of (fT7|) . Equivalently, we need to show that 

To verify this, write £7 = {iti < • • • < u m }, S 10 = {s^ < • • • < s' n _ m }, and let X be the ordered 
sequence (u±, . . . , u m , s[, . . . , s^_ m ). Consider the determinant 

D= det «*,/,-». 

«£A,l<j<n 

On the one hand, reordering the rows gives 
§m(m+l)+ £ s 

D = (-l) ^ fe/ det {{eijj)). 

l<i,j<n 

On the other hand, applying Laplace expansion to the first m rows gives 

Eim(m+1)+ J2 t 
(-1) ^ T det (lei,fi)) det (( ei , /,-)). 

TC[n],\T\=m 

This completes the proof. ■ 
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Consider the case of Lemma [3] when &i(x) = fi(x) = x % . Then, es{x) = A(x)s\(x), where 
S = {si < ■ • ■ < s m } and Aj = s m +i_j — i. Thus, noting also that 

n 

det {(etjj)) = det = TT(pi_i,pj_i), 

we obtain Proposition [5J 
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